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ABSTRACT 

A 13°xl3° image of Balmer Ha emission at galactic latitude -65° is presented. 
Sensitivity is limited in part by confusionll: the peak-to- valley anisotropy of the Ha 
surface brightness is ~0.2 Rayleighs on angular scales of 0.1° to 1.0°. The morphology 
of the Ha emission is similar to that of the 100 /jm emission previously observed by 
the Infrared Astronomical Satellite (IRAS). A point-by-point comparison of the two 
shows a marginally detected (3a) positive correlation, p = +0.14±0.04. The slope of 
the correlated components of the 100 fim and Ha emissions is 1.26]to32 MJy sr _1 R _1 . 
Using parameters from the literature, we predict that emission from dust and ionized 
gas at high latitudes produces ~3 times more Ha emission per unit 100/im emission 
than does backscattering of Ha emission from Galactic H II regions by dust at high 
latitude. Observations of this type may allow us to distinguish between Galactic 
foreground and cosmic background for both the free-free emission and the thermal 
dust emission associated with the warm ionized medium of the Milky Way. 



Subject headings: 

1. Introduction 

Martin Harwit ( Harwit 1981| ) might say we have rediscovered in Ha the interstellar cirrus 



that was discovered first in the infrared by IRAS ( Lo et al. 1984| ). Such a rediscovery is not 
insignificant. 

To measure anisotropy of the cosmic microwave background, astronomers must first measure 
and subtract foreground emission from the Milky Way's interstellar medium (ISM). Three 



1 Alfred P. Sloan Research Fellow. 



The Ha surface brightness at high Galactic latitudes 
1 R = 10 6 /4tt photons cm" 2 s" 1 sr" 1 ( |Reynolds 199^ ). At Ha 1 R = 
measure EM = 2.75 cm _6 pc for gas at 10 4 K. 



is approximately 
= 2.42 x 10~ 7 ergs cm" 



1 Rayleigh (R), where 



or an emission 
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components of foreground emission from the ISM are important (Bennett et al. 1992): 1) 
bremsstrahlung (or "free- free") emission from free electrons accelerated by interactions with ions 
in a plasma, 2) synchrotron emission from relativistic electrons accelerated by magnetic fields, and 
3) thermal emission from dust. The Ha images directly trace the free-free component, because 
ionized gas emits Ha and bremsstrahlung in direct proportion. In addition, the Ha images 
indirectly trace the thermal emission from that dust associated with the ionized medium. 

This paper shows that Ha and 100 fim emissions are weakly correlated at high galactic 
latitude. High latitude, filamentary nebulosities have been observed before ( |Guhathakurta fc Cutr 



1994| , pandage 1976| , and references in each), but verifying the emission to be Ha in particular had 
not been possible. In fact, Guhathakurta &: Tyson (1989) showed that the optically broad-band 
red color of cirrus is not due to Ha emission, but may be due to luminescence of hydrocarbons. 
Kogut et al (1996) have shown that multi- frequency radio emissions (commensurate with free- free) 
are correlated with 140 fim emission. Broadbent, Haslam, & Osborne (1989) have shown that 
free-free emission and thermal dust emission are positively correlated at low galactic latitude, i.e. 
emissions from H II regions. 

At low galactic latitudes, there is considerable confusion for optical tracers of gas: where 
there is gas, there tends to be dust too, and the dust that is emitting at 100 //m, for example, is 



absorbing Ha light. Even at high latitudes, where dust opacity is minor ( Reynolds 1980 ), for any 
integral along a line of sight, the 100 (im. surface brightness is proportional to the density of dust, 
whereas the Ha surface brightness is proportional to the electron density squared (see Section 
H). For these reasons, we do not expect a strict proportionality between the Ha and infrared 
emissions. 

Our Ha images provide not absolute surface brightnesses but only relative ones. The same is 
true of the IRAS images. We do not attempt absolute calibration of the Ha surface brightness 
because our Ha filter's 1 nm bandpass is too wide to suppress geocoronal Ha, which is typically a 
few Rayleighs and varies with time (Broadfoot & Kendall 1968). Fabry-Perot spectroscopy allows 



absolute calibration of Ha surface brightness, and the infrared has been calibrated by COBE. 
However, to compare IRAS and Ha images, it suffices to filter out the background in the same 
way for both images, and that we have done (see Section ||). 



2. Instrumentation and Observations 

The observations were made at the Mt. Laguna Observatory with a wide-field, narrowband 



imaging camera similar and yet superior to other such cameras (e.g. Bothun fc Thompson 1988] ) . 
We attached it "piggy-back" upon a 1-m telescope, which was used solely for blind pointing 
and tracking; no acquisition or guiding was necessary. The camera consists of a filter wheel, a 
135-mm f/2.8 lens, and a 2048x2048 pixel CCD. With the 135-mm lens, the 15-/im pixel spacing 
corresponds to 0.38'/pixel, and a field of view of 13°xl3°. 
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To tune the center wavelength of the 1.0 nm wide (FWHM) Ha filter, the filter wheel was 
thermoelectrically cooled and thermostatically controlled to a precision of 1° C. The temperature 
was set such that the Ha filter's transmission of zero-velocity Ha peaked off the optical axis while 
maintaining the on-axis transmission to greater than 75% of the peak. At off-axis angles of > 4.8°, 
the transmission of zero- velocity Ha drops below 50%, so each image is cropped at a radius of 4.8°. 
Other than cropping, we made no attempt to correct for the radial dependence of the sensitivity 
to zero- velocity Ha .Considering that and other potential errors, we estimate our absolute scaling 
of Ha surface brightness could be in error by as much as a factor of 1.5. Observations of external 
galaxies show the velocity coverage extends from —250 km/s to +250 km/s. 



The backside illuminated CCD has a quantum efficiency at 656.3 nm of 95% ( Leach 1996| ) 



The sensor is cooled with liquid nitrogen to -120° C, so for these observations, dark current is 
negligible. The noise of the readout amplifier is 9 e~ rms, but the integration times were selected 
so that the readout noise would be negligible compared to the Poisson noise of the sky's photons. 

Forty-eight pairs of exposures were taken over a period of 5 nights, September 21-25, 1995. 
Each pair consisted of a continuum image followed by an Ha image while the telescope turned 
at the sidereal rate for the duration of the pair. Each continuum image was a 1-minute exposure 
through a R band filter, except for the first few, which were 2 minutes each (and which had too 
many saturated stars). Each Ha exposure was 20 minutes; the total observing time on-source 
was 16 hours. The images were dithered; that is, the telescope was repositioned after each pair of 
images. The dithering pattern was a 5x5 grid with 1° spacing, covered almost twice with the 48 
image pairs. 

We observed the North America nebula and Barnard's Loop to calibrated the Ha surface 
brightness. These two H II regions have surface brightnesses of 850 R and 153 R (each ±10%, 



Reynolds 1980 ), averaged over 48' diameter fields centered on 21h00m.0 +45°00' [2000] and 



05h53m.8 -06°08' [2000]. 



3. Data Reduction 

Much of the data reduction is aimed at eliminating stars from the images while preserving 
the diffuse Ha emission. The stars are suppressed in a number of ways: 1) continuum subtraction, 
2) flagging pixels that are bright in continuum, 3) flagging pixels at the locations of filter-induced 
flares near very bright stars, and 4) median filtering with a spatial scale larger than a stellar image. 

The CCD bias was determined by overscanning. The relative gain correction for the CCD 
pixels (the "flat field") was determined by observing the zenith at sunset in order to minimize 



gradients and curvatures ( Chromey &: Hasselbacher 1996 ). Each continuum image was scaled by 



a constant and subtracted from the Ha images. The constant was empirically-determined based 
upon stars in each field. 
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The continuum-subtracted images invariably exhibited a smooth intensity pattern with radial 
symmetry about the optical axis. The radial pattern may be due to telluric line emission, although 
it could also be due to the combination of vignetting and inaccuracy of the flat fields. In any 
event, the circular pattern was removed by fitting and subtracting circular isophotes. 

Our experience has been that continuum subtraction is adequate for faint stars, but 
inadequate for bright stars. Invariably, the two images are not focused identically or they are 
not registered well enough or both. In the presence of these systematic errors and the increased 
Poisson noise at the positions of bright stars, it is worthwhile to excise those pixels near bright 
stars from the data. We flag pixels in the Ha images based upon their respective intensities in the 
R band images, and flagged pixels are excluded when computing the medians described below. 

Each image was median filtered spatially with a 2° diameter ring median filter ( pecker 1995| ). 
Each median-filtered image was subtracted from its original to eliminate large scale intensity 
variation due to the variation in airmass across the field, residual inaccuracies in flat fielding, or 
the occasional telluric cirrus cloud. To check that this procedure was not eliminating real emission 
from the Milky Way, we also reduced the data without this median filtering step, but instead 
substituted fitting and subtracting a plane of intensity. The two methods give similar results. 

The images were registered automatically to a common coordinate system, using a few dozen 
SAO stars on each image to determine the appropriate shift, expansion, rotation, and distortion. 
The 48 registered images were combined by averaging them pixel by pixel, while excluding flagged 
pixels. Also, an iterative procedure eliminated those individual pixel values that were further 
than 2a from the median of the other values for that particular pixel. Finally, a ring median filter 
of radius 2.3' was applied, and then the image was convolved with a circular Gaussian of width 
a = 1.5'. 



4. Results 

The Ha image (Figure |l|) appears similar to the corresponding area of sky observed by IRAS 
at 100 iim. (Figure |2|). We have inscribed boxes in Figures || and || that contain the two most 
apparent similarities: 1) at left, a filled arc of 100 /im emission, enveloped by an arc of Ha 
emission, and 2) at center-top, a potato-shaped outline in both 100 /im and Ha emission!! 

To quantify the association between the Ha and 100 /im emissions, we evaluated their 
(Pearson) correlation coefficient, p = +0.115, over the central 8.5° x8. 5° subregions of the 13°xl3° 
images. The sugregion was selected because the effective observing time is less near all edges of 
the Ha mosaic and because residual, circular artifacts of the mosaicing process are apparent near 
the western edge. 



' A "movie" that "blinks" the two Ha and 100 /im images is available from the author. 
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We also correlated the Ha image with 36 controls. The controls are 100 pm IRAS images 
from "random" but statistically similar directions on the sky, centered at galactic coordinates b = 
±67°, 1 = 1°, 21°, 341°. Each control was created and processed like the image in Figure 2. The 
controls should be uncorrelated with the Ha image (i.e. the expectation value, E(p) = 0), so they 
provide a standard by which to judge the validity of our measured correlation. The distribution 
of correlation coefficients of the 36 controls is approximately Gaussian, with p = -0.017 and a = 
0.041 (see Figure |). 

The measured correlation coefficient (p = +0.115) is statistically significant from zero by 
+2.8(7. From the controls we can conclude with ~95% confidence^ that p ^ 0, but we should not 
immediately conclude that p = +0.115+0.041. Mathematically, we expect p to be larger than 
+0.115 because noise in the data tends to reduce \p\ from its true value. Physically, we expect p 
to be larger or smaller for various directions on the sky. For example, in the boxes of Figures 1 
and 2, which were inscribed because the emissions within them looked similar, p = +0.150. 

The bias of p due to noise in the two images can be estimated, and it is not large. We let of 
equal the variance of the infrared intensity and we let a\ equal the variance of the Ha intensity, 
both as measured over the inner 8.5° x 8. 5° subregions. The af and o\ each include both genuine 
variation of emission and random noise. We evaluate the noise variances in each image, r/f and 
r]h, by computing the average of the variances in small regions apparently lacking emission. If 
the ratios (xf/r/f and cr^/rj^ are both large compared to unity, the bias is small and more precise 
estimates of the noise are not warranted. The measured values are a\lr\\ = 2.4 and UhMh = 2-6. 
In can be shown that the ratio of the true to the observed (Pearson) correlation coefficients is 



which in this case is 1.19 and is smaller than the fractional uncertainties of both the statistical 
variation of p estimated from the controls and the spatial variation evident from the value 
of p = +0.150 of the alternative subregion. Accounting for the bias of noise, we estimate 
p = +0.14+0.04. 

Having established that the 100 pm and Ha emissions are weakly correlated, we evaluate 
the ratio of the two correlated emissions by a linear fit of the form, 1(100 pm) = a I(Ha) + b, to 
the data within the central 8.5° x8. 5° subregions of the two images. The slope "a" is of interest; 
the value "b" is not, because the zero points of both the 100 pm intensity and the Ha intensity 
are arbitrary. Because both the 100 pm and Ha emissions have substantial observational errors, 
ordinary least squares fitting is invalid .i We fit the data by two methods as described below. Both 



4 That 2 of 36 controls share the same bin in Figure 3 as the true correlation coefficient indicates a ~5% likelihood 
that the two emissions are uncorrelated (p = 0), i.e. that the measured correlation is due solely to chance. 

5 Inthis context "ordinary" least squares refers to fitting a "dependent" (noisy) variable Y against an "independent" 
(noise-free) variable X. 




(1) 
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methods allow for errors in each of the observed quantities. 



The first method is based upon the work of Feigelson & Babu 1992 B It yields a value for 
the slope a = 1.2(±0.12) MJy sr _1 R _1 , under the assumption of additive Gaussian errors of 
zero mean and variances of (0.17 MJy sr -1 ) 2 and (0.077 R) 2 , for the 100 fim and Ha emissions, 
respectively. We determined the variance of the Ha emission, (0.077 R) 2 , from an ensemble of 
1-dimensional Gaussian fits to the bivariate distribution of 100 /im and Ha emissions by rows; 
the variance of the 100 /xm emission was determined similarly, but by columns (cf. Figure 4). 
Note that those variances include real (but uncorrelated) emission in addition to detector and 
shot noise; the uncorrelated emissions act as approximately Gaussian "errors" with respect to the 
correlated components of the emissions. We consider the formal error (cr a = 0.12MJy sr" 1 R" 1 ) to 
be a lower limit of the true error. To investigate empirically the reliability of the slope estimate, we 
performed the same analysis separately on the four quadrants of the central 8.5° x 8.5° subregion. 
The slopes so measured were 1.5, 1.4, 1.3, and 0.6 MJy sr" 1 R _1 for the NW, NE, SW, and SE 
quadrants, respectively. 



The second method is a simplification of the "BCES" method of Akritas &: Bershady 



1996. BCES refers to "bivariate correlated errors and intrinsic scatter"; our simplification is 
to assume that the observational errors of the 100 fim and Ha emissions are uncorrelated. We 
evaluate the observational errors by examining small regions apparently lacking emission. The 
standard deviations so measured are 0.07 MJy sr" 1 and 0.04 R, for the 100 /jm and Ha emissions 
respectively. The BCES bisector slope (see Akritas fc Bershady 1996| ), evaluated over the central 



8.5° x8. 5° subregion, yields the same slope as the first fitting method, a = 1.2 MJy sr _1 R _1 . 
Likewise, the BCES bisector slopes for the NW, NE, SW, and SE quadrants are 1.68, 1.31, 1.16, 
and 1.53 MJy sr _1 R , respectively. 

The result of averaging the 8 slopes from the four quadrants and the two methods 
is a = 1. 26^Q 32 MJy sr _1 R _1 . The mean value given is a vector average (i.e. 
< a > = tan(< arctan(a;) >), i = 1, 8). The excursions (+0.45 and —0.32) describe 
the la dispersion of the 8 slopes, also calculated vectorially, not arithmetically. Additional 
observations are required to determine if the dispersion of slopes is statistical or intrinsic to the 
ISM. 

The brightest object in the Ha image has no counterpart at 100 ^m, illustrating that the 
correlation between 100 fim and Ha emissions is not true in general. To the contrary, additional 
observations at lower galactic latitudes (|6| < 30°) and of less sensitivity show many Ha features 
with no infrared counterparts ( McCullough 19971 ) . Those additional images also show Ha-infrared 



correlations, both positive and negative. At the lower Galactic latitudes, there is considerable 
confusion. 

The brightest object in Figure || may be associated with the Magellanic Stream; it lies 



'Note that their variable I — dy/a^, not u^/^y as originally printed. 
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between two H I peaks of MS IV ( Cohen 1982| ) and its long axis is parallel to the Magellanic 



Stream. If so, it is the brightest Ha emission from the Magellanic Stream (cf. Weiner & Williams 
1996, Weiner 1996] ) . However, it might be simply a density enhancement in the ISM unrelated to 



the Magellanic Stream. The critical test is to measure its velocity. 



5. Discussion 

The Ha light we observe could be from either a reflection nebula or an emission nebula, or 
both. In this Section we show that the reflection component is expected to be small compared to 
the reflection component. 

In the case of a reflection nebula, the surface brightness of the Ha light is proportional to 
the product of the surface emissivity in Ha of the disk of the Milky Way and the dust column 
density along the line of sight. Using canonical values for the surface emissivity, the ratio of dust 
to gas, and the dust's scattering cross section and phase function, Jura (1979) predicts the surface 
brightness of the Ha light I Q (in Rayleighs) backscattered from dust: 

N H 

^ = °- l0 1n20 -2 > ( 2 ) 

l(r u cm £ 

where Nh is the total hydrogen column density. 

In the case of an emission nebula, the surface brightness of the Ha light I a (in Rayleighs) is 
proportional to the emission measure QReynolds jggg ): 



I Q = O.36T4 a92 EM, (3) 
where T4 is the electron temperature in units of 10 4 K. In the warm ionized medium of the 



ISM, T4 ~ 0.8 ( Reynolds 1992 ). The emission measure EM = J n e n p dl ~ / n^dl has units 
cm -6 — pc. The region is assumed to be optically thick to Lyman lines (a.k.a. Case B), and 
extinction is ignored, but otherwise, equation || is derived directly from the fundamental physics 
of recombination. 

From the emission measures and dispersion measures along the high-latitude lines of sight 
to a few globular clusters containing pulsars, Reynolds (1991a) shows that the density of high 



latitude ionized gas, where it existsfj is ~0.08 cm" 3 and the ratio of the ionized hydrogen column 
density to the total hydrogen column density ranges from 0.2 to 0.4 ( [Reynolds 1991b] ), with an 
average of 0.27 for the five lines of sight. 

Combining the emission and reflection contributions to the Ha surface brightness using the 
above values from the literature, we find 

-0.92 N H+ , n 1 n\ N H 



I a = (11.7T4 U - W ne^-+0.10) 



4 e N H ; 10 20 cm- 2: 



7 Reynolds estimates that the ionized gas fills ~20% of the volume. 
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= (°- 31+0 - 10 )l020^' ^ 

where we have expressed the EM as a product of the volume density, n e , and the ionized hydrogen 
column density, N#+, by approximating J n e n p dl ~ n e J njr+dl. The latter approximation is valid 
if hydrogen is the dominant ion and if the electron volume density is bimodal, i.e. either n e is 0.08 
cm -3 or it is zero. 

Equation ^ predicts the emission component to be ~3 times greater than the backscattered 
component. The backscattering component will be less than a third of the emission component 
if the gas is either denser or more highly ionized than we have assumed. Reynolds et al (1973) 
estimate the ratio of backscattered to total Ha light to be 0.1 (—0.1, +0.2) toward a few Lynds 
nebulae. 



It is well known that dust and neutral hydrogen are positively correlated in general ( [Boulanger 



fc Perault 1988j ), but is that true in this region? The three H I maps of highest resolution and 



sensitivity ( |Hartmann 1994 ,Heiles & Habing 1974, Cohen 1982] ) are inadequate to allow us to 



answer this question because the beams are too large (0.5° FWHM) and the sensitivity is too 
poor (the r.m.s. is equivalent to Nh ~ 1 to a few times 10 19 cm -2 ). The peak 100 /mi emission, 
~0.5 MJy sr -1 above the background, corresponds to a peak H I column density of 6xl0 19 



cm -2 , assuming Iioo = 0.85NH/10 20 cm~ 2 MJy sr _1 ( Boulanger fc Perault 1988 ). Nevertheless, we 
analyzed the H I column density data of Heiles & Habing (1974) exactly as we did the Ha data 
above and found no statistically significant correlation with the 100 /mi data.i The H I&infrared 
correlation coefficient is —0.119, which corresponds to 1.5a less than zero (a = 0.078 is evaluated 
as before using thirty six 100 /im controls). The anticorrelation of H I&infrared in the same region 
where Ha&infrared is positively correlated is intriguing, but it is not statistically significant. The 
anticorrelation also is opposite in sign to the firmly established results of Boulanger & Perault 
(1988). 

In Section ^ the slope of the correlated components of the 100 /mi and Ha emissions is 
estimated to be 1.26±g:f| MJy sr" 1 FT 1 . From 4 years of COBE data, Kogut et al. (1996) 
determine the slope of the correlated components of the antenna temperature Ta of the 
DMR at 31.5 GHz with respect to the 140 fim surface brightness I140 measured with DIRBE: 
T A /Ii40 = 6.37(±1.52) /iK(MJy sr" 1 )" 1 . In so far as free-free emission dominates thermal 
dust emission at 31.5 GHz, the slope 6.37(±1.52) /iK (MJy sr -1 ) -1 correspondsi to a ratio 
Ii4o/Ihc* = 0-90lo; 29 MJy sr _1 R _1 , which is similar to the Iioo/Ih« ratio derived in this work. 



On the ~7° angular scale of DMR, Kogut et al. 1996 conclude that more than one third of (and 



perhaps "the bulk of") the free- free emission must be from that component correlated with 140 
/tm emission. At much smaller angular scales, the bulk of the Ha emission is not correlated with 



8 The data of Cohen 1982 



exist only on 8-hole paper tape, and the data of Hartmann 1994 were not yet publicly 



available in machine readable form when this paper was written. 



See e.g. Bennett et al. 1992 equation 4. 
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100 /xm emission (Section ||); thus it appears that the correlation between ionized gas and dust 
emissions increases with increasing angular scale. 

6. Summary 

An Ha image of a 13°xl3° region at (l,b) = (71°,— 67°) shows structure ~10 times fainter 
than the 12°xl0° region of sky at (l,b) = (144°, -21°) previously observed by Reynolds (1980). 
The new image gives a preliminary view of the sky over many degrees at unprecedented sensitivity. 
The new Ha image weakly resembles the 100 /im IRAS image; a point-by-point comparison of the 
two shows a positive correlation (p = +0.14 ± 0.04); and comparison with COBE data indicates 
that the correlation coefficient increases with increasing angular scale. The slope of the correlated 
components of the 100 ^m and Ha emissions is 1.26^g' 32 MJy sr _1 R _1 . 

An important extension of this work (already underway by many investigators) will be 
to image the Ha emission over much of the sky. An all-sky Ha image may prove critical to 
distinguishing the Galactic foreground from the cosmic background at the precision necessary to 
determine cosmological parameters. 

Jeremy Hicks assisted with the observing. Bob Leach made the CCD camera. We appreciate 
discussions with John Gaustad, Carl Heiles, Margaret Meixner, William Reach, Ron Reynolds, 
and Dave Van Buren. We thank Ben Weiner for providing data in advance of publication. 
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Fig. 1. — This Ha image is 13° on a side and shows features as faint as 0.2 Rayleighs. The triangular 
object near the center is the brightest object in the field, at 0.6 R. The inscribed box contains two 
Ha emission features associated with IRAS 100/um emission: the arc at left and the potato-shaped 
outline at top-center (cf. Figure 2). The center of the image is 23h40m00s -12 o 00'00"[2000], (l,b) 
= (71°. 5, —67°. 4). North is up and East is to the left. The Galactic plane is toward the upper-right 
(NW). 
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Fig. 2. — This IRAS 100/im image shows the same region as Figure 1. Because of its proximity to 
the ecliptic, this field is one of the so-called "reject" fields of IRAS. To reduce variations on large 
angular scales, a 2° diameter ring-median filtered version of this image was subtracted from the 
original. The same filtering was applied to Figure 1. 
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Fig. 3. — This histogram of the absolute values of the correlation coefficients of 36 controls shows 
that the correlation coefficient of the Ha and WOfim emissions in Figures 1 and 2, +0.115 (arrow), 
is statistically significant at 3a. 
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Fig. 4. — The 100/xm versus Ha bivariate distribution is plotted (contours); the best fit line (solid) 
and lines with ±la deviations of the slope (dashed) are superimposed. The contours are linearly 
spaced at intervals of 12.5%, 25%, ... of the peak; the 50% contour is bold. The data are from the 
central 8.5° x8. 5° subregions of the 13°xl3° images in Figures 1 and 2. 



